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Abstract: Flic main characteristics oi the Oner, us pollination season were studied in 14 
different localities of the Iberian Peninsula from 1992 2004. Results show that Quercus 
flowering season has tended to start earlier in recent years, probably due to the 
increased temperatures in the pre-flowering period, detected ai study sues over the 
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developed by the Hadley Meteorological Centre, in order to determine the expected 
advance ni ihe start of Qiwrctn pollination in future years. At each stud) site, airborne 
[iolle.ii curves picsented a similar pattern in all study years, with different peaks over the 
season attributable in many cases to the presence of several species. High pollen 
concentrations were recorded, particularly at Mediterranean sites. This study also 
proposes forecasting models to predict both daily pollen values and annual pollen 
emission. All models were externally validated using data for 2001 and 2004. with 
acceptable results. Finally, the impact of the highly-likely climate change on Iberian 
Quercus pollen concentration values was studied by apply m:' R( AI meteorological data 
lor different future years. 2025. 2050, 2075 and 2099. Results indicate that under a 
doubled C0 2 scenario at the end of the 21st century Quercus pollination season could 
start on average one month earlier and airborne pollen concentrations will increase In 
50'/r with respect to current lev els. w nh higher values in Mediterranean inland areas. 
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INTRODUCTION 

In spring, tree pollen grains of the Fagales order which 
includes Quercus genus, are a mam allergen source in the 



northern hemisphere, with cross-reactivity oeeunin 
between pollen within the same order [40. 44]. Severn 
studies consider Quercus pollen type responsible lot 
allergy in areas with abundant Quercus vegetation |5. 38 
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46, 52, 56, 57]. The Que reus genus is represented in the 
Iberian Peninsula by 1 1 tree species commonly known as 
oaks (cork-oaks, holm-oaks, gall-oaks. etc.). These 
species are distributed throughout the entire peninsula, are- 
well represented in Atlantic forests and dominate most 
Mediterranean ecosystems. All are anemophilous trees 
producing stenopalynus pollen grains mot distinguishable 
under light microscopy) in high quantities that are well 
dispersed through the air 129. 58 1. The high abundance of 
these species in natural and semi-natural vegetation areas 
increases allergy risk for agriculture and forestry -related 
workers in the Iberian Peninsula, highly exposed to 
Quereus pollen grains [56], Que reus pollen season 
dynamics and their relationship with meteorological 
factors were studied at 14 sites of the Iberian Peninsula 
(Spain). Particular attention was paid to the effect of the 
detected temperature increase in later years on the 
Quereus pollen season, and to the potential consequences 
of the global climate change. 

In recent years, the rise in winter and spring tempera- 
tures as a consequence of global climate change is 
causing earlier leafing and flowering in many species, in 
both Europe and the USA [4, 17. 41]. Changes in the 
pollination period features, both tuning and intensity, of 
plants producing allergenic pollen would have important 
consequences on allergy sufferers [3. 16]. Knowledge and 
forecasting of the pollination dynamics provides useful 
information for managing and preventing allergic symp- 
toms. The effect of recent climate warming on vegetation 
phenology and aerobiology is an important issue which is 
being taken into account in climate-change studies [28, 
53]. Because Quereus pollen emission is a springtime 
phenomenon, pollen season dynamics are sensitive to 
changes in phenology. 

Global average surface temperature has increased ().(>"(' 
since the late 19th century. Spring temperatures in the 
Mediterranean countries have risen by 2°C over the last 
two decades [36, 42]. The increase has been shown to 
prompt earlier leaf bud bursting and flowering in 
temperate Mediterranean tree species such as oak and 
olive [24. 27 1. The Special Report on ['missions Scenarios 
(SRES) produced a series of scenarios of C0 2 emissions 
representing outcomes of distinct narratives of economic 
development [45 J. Scenarios arc plausible future states of 
the climate system. They are derived from predictions of 
climate change from Global Climate Models (GCMs) 
based on CO; emissions scenarios. GCMs are mathema- 
tical descriptions of the phy sical elements and proc esses 
in the atmosphere, oceans and land surface (winds and 
ocean currents, clouds, rainfall, soils). However, in areas 
where bio-geography has a significant effect on wheather, 
scenarios based on global models will fail to capture local 
details for impact assessment. The best method for adding 
these details to global models is to use a Regional Climate 
Model (RCM), essentially a higher-resolution version of a 
GCM covering a limited area [48]. Simulation RCM can 
be replicated by using forecasting phonological models, 
and allows assessment of the potential variation in 



vegetation dynamics [47]. Here, we examine the potential 
changes on Quereus floral phenology in Spain, indicated 
by airborne pollen dynamics. We use the findings of 
recent studies and the models developed in this study to 
drive them with replicated RCM simulations, in order to 
explore potential changes in Quereus pollen season that 
could accompany future temperature changes in the 
Iberian Peninsula. 

MATERIAL AND METHODS 

Description of study areas. The 14 study sites present 
different climatic, topographic and vegetation characteris- 
tics: some of these are summarised in Table 1. Northwest 
(N'W) Spain is a mild-temperature, rainy area influenced 
by the Atlantic Ocean. The predominant form of vegeta- 
tion is deciduous woodland, with a large population of 
Quereus species. Vigo and Santiago are located in the 
Euro-Siberian region, with similar weather and vegetation 
conditions and low altitude. Quereus robur L. (first in 
flowering) and Quereus pyrenaiea Willd. are the most 
common Quereus species in the forests of this area. The 
remaining study areas are situated in the Mediterranean 
Region, although Ourense and Leon are highly influenced 
by the Euro-Siberian region climate. The sampling sites 
located in the northeast (NE) of the Iberian Peninsula, 
Barcelona, Tarragona, Lleida and Girona, have a typically 
Mediterranean landscape, with, extensive areas of holm 
oaks (Quereus ilex subsp. ilex L., Quereus ilex subsp. 
ballota (Desf.) Samp), Quereus humilis Miller, and 
Quereus suber L; Kermes oak (Quereus coccifera L.) 
woodland appears on the most impoverished soils. The 
climate of the first 2 sites is characterised by the 
proximity of the Mediterranean sea. Temperatures are 
warm and annual rainfall is scarce. The phenologica) 
flowering order is as follows: Q. humilis, Q. coccifera, Q. 
suber and Q. ilex subsp. ilex and Q. ilex subsp. ballota 
The Central area has a dry, continental climate with 
significant inter-seasonal variations in both temperature 
and rainfall. In the local countryside around Madrid, the 
dominant Quereus species is the holm-oak Q. ilex subsp 
ballota, which forms large areas of oak-grass savannahs 
used for livestock farming and known as dehesas. Leon, 
located to the north, is influenced by the Euro-Siberian 
region. Its landscape is characterised by mid-altitude holm 
oak woodland, although higher woodland populated by Q. 
pyrenaiea and Q. faginea Lam. can also be found. These 
2 species flower ) month later than the holm-oak. The 
Southern region includes 5 of the study sites: Cordoba. 
Priego. Jaen, Granada and Malaga This region has a 
warmer, drier Mediterranean climate, with more moderate 
conditions in the coastal city of Malaga. Typical 
vegetation includes dehesas of holm-oak ( first flowering), 
cork-oaks (last flowering) in the most humid areas with 
acid soils, and kermes-oak (second flowering) in 
impoverished soil areas. A few populations of Q. 
pyrenaiea and Q. faginea grow in deep ravines and humid 
locations in the Supramediterranean vegetation area. 



etting in the Peninsula, elevation (me 
cJ data and Oucrcus species present at 

Coordinates Mean air temp. W. 



North-West Santiago 

North-West Vigo 

North-West Ourense 

North-East Barcelona 

North-East Tarragona 



7 years Q. robur (L.) 

i 1 995 200 1 ) Q. pxrenaica (Willd.) 
O. suberlL.) 

9 years Q. robur (L.) 
( 1 993-200 1 ) Q. pxremuca (Willd. ) 
Q. sufcei" (L.) 

Q. ilex subsp. bailout (Desf.) Samp. 

8 years Q. ilex subsp. //« (L.) 
(1994-2001) £>. (7c.v subsp. /.„//<»(<( (Dest.) Samp. 

Q. humilis (Miller) 
(2. cocci/fern (L.) 
Q..vufee;(L.) 

0. cerrioides (Willk. & Costa) 
f> years (_>. i/c.v subsp. i/ci (I..) 
(1996-2001) 0. i/c.v subsp. f> ( i//o/<i (Desl'.i Samp. 
Q. humilis (Miller) 
(A pyrenaica (Willd.) 
Q. faginea (Lamk) 



North-East Girona 



( 199(,-2001 I O ilex subsp. 

Q. humilis (Miller) 



North-East Lleida 
Centre Leon 



. ballota (Desf.) Samp. 



9 years (J. ilex subsp. ballota (Desf.) Samp. 
(1993 2001 ) Q. prreiuaca (Willd) 

Q. ronrifrw (L.) 
Q.faginea (Lamk) 

10 years Q. t/e.v subsp. fc«//ofa (Desf.) Samp. 



South 
South 



3 years Q. ilex 
-2001) (2>. «><• 
> Q.ilex 



era 



ballota ( Desf ) Samp. 



t subsp. ballota (Desf.! Samp 
(1996-2001 Q. coreifera (L.) 

Q. faginea (Lamk) 
Id scars (J. iie\ subsp. ballota iIVsl' i Samp 
(1992-2001) 0. eoeeiferet (L.) 

Q. .vh/kt (L.) 



1 0 yea 



v subsp. />«//()(» (Desf.) Samp 
' i/Wv, (L.) 



Airborne data. The aerobiological behaviour of 

Quercus pollen was statistically analysed using a 13-year 
database (1992-2004) of pollen and meteorological data 
from 14 sites in the Iberian Peninsula (Tab. 1, Fig. 1). 
These data were drawn from the Spanish Aerobiology 
Network (REA) data bank held at the REA Coordinating 
Centre (University of Cordoba). In all the stations. Hirst 
volumetric spore traps were used. They are located on 



buildings approximately 15 m above ground level. These 
traps have an autonomy of 7 clays and collect pollen 
continuously with a given absorption flux, enabling daily 
and even hourly pollen concentration data to be obtained. 
Daily mean data are expressed in terms of pollen 
concentrations per cubic metre of air (gr/nr). The 
aerobiological methodology recommended by the REA 
was used [ 15]. Pattern curves were constructed w ith 5-day 
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Figure 1. Predicted temperature increase lor 2()<« m the Iberian 
Peninsula (Regional Climate Model HadCM3A2a model). The study 
localities are indicated on the map 

moving-average values calculated from daily values. 
Graphs were plotted following the line-up method [32J. 
The Pollen Index (PI) was defined as the total pollen 
detected in the Main Pollen Season (MPS) of a given 
year. The start of the MPS was defined as the clay on 
which 1 pollen grain/m' was reached wherever 5 sub- 
sequent days contained 1 or more pollen grains/nv'. The 
end was determined as the last day on which 1 pollen 
grain/m 3 was recorded wherever subsequent 5 days 
presented concentrations below this level [25]. Data from 
1992-2000 were used to build prediction models. These 
were validated in 2001 and 2004. except in L.eon and 
Priego. In Leon, models were validated in 2000 (not 
included in the model for Leon) and 2004. due to 
technical problems that interrupted sampling during the 
spring of 2001. In Priego, models were validated in 2001 
(not included in the model for Priego) and 2002. 

Meteorological data. The Spanish National Institute of 

Meteorology (wuu.inm.es) supplied ilaih meteorological 
data from its stations located nearest to the pollen traps. 
The Iberian Peninsula suffered a severe drought between 
1990-1995, mainly affecting Mediterranean \treas: as a 
result, the data were divided in 2 groups for the statistical 
analysis (except for the NW sites, unaffected by the 
drought). The extremely dry period. 1992-1995. and the 
usual climate period. 1996-2000. 

Statistical analyses. The STATISTICA© software 
package was used in till analyses. For forecasting daily 
and annual pollen concentrations, a Spearman correlation 
test was performed previously to detect the chief variables 
affecting pollen values. Stepwise multiple regression 
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analysis was performed in order to construct both daily 
and annual pollen forecasting models. Predictive variables 
were chosen by applying a Factor Analysis for Principal 
Components in order to avoid multi-colinearity. For daily 
forecasts, previous attempts b\ analysing the whole- 
season revealed very weak and contradictory results, 
therefore 2 types of regression models were developed: 
Pre-Peak models Groin the start date up to the maximum- 
value date): and After-Peak models (from the maximum- 
value date up to the end of the pollen season). Daily 
models were developed by constructing general equations 
comprising several years. For Mediterranean sites, 2 
different groups of years were taken into account, to 
reflect the drought period (1992-1995; 1996-2000). Daily 
variables, such as weather data from the previous day and 
accumulated weather data from the previous 2-5 days, and 
previous and cumulative pollen counts up to the date for 
which the forecast was made were used as predictive 
variables. For PI forecasting models, the following were 
calculated: fortnightly minimum, maximum and mean 
temperature; fortnightly total rainfall; fortnightly evapo- 
transpiration and chilling hours based on daily data, to be 
used as independent variables in order to obtain predictive 

Future start dates were calculated by applying future 
meteorological data to a previously-developed thermal 
forecasting Growing Degree Days (GDD") method [27]. 

All predictive models were validated in years not 
included in the equations. For the daily model, 2 different 
equations were developed for Pre-Peak and After-Peak 
pollen curves. Since weather conditions from 2000-2004 
were similar to those of 1 996-2000 in the Mediterranean 
area, the proposed forecast formula for that period was 
used for validation purposes at years 2001 and 2004. 
Expected vs. actual data were compared graphically. The 
R2 coefficient and the Wilcoxon matched pairs,, test were 
performed between expected and actual data. Annual 
pollen models results were validated by calculating Root 
Mean Square Error and % Error. 

Simulation model. Future meteorological data model- 
led by the Hadley Meteorological Centre including in the 
Intergovernmental Panel on Climate Change (1PCC), 
were used for estimating future pollination start dates. 
The data were obtained through a contract (Contract EPG 
1/1/124) signed by Cordoba University and the Climate 
Impacts LINK Project (Hadley Centre). The present study 
also analyses the possible future consequences of this 
expected climate change on the annual amounts of 
released Que reus pollen, based on the forecast equations 
developed here. The model is a mathematical description 
of major elements and processes in the atmosphere, 
oceans and kind surface comprising the climate system. 
By applying statistical downscaling the method adds fine- 
scale information. This technique uses observation in 
historical meteorological data sets and today's climate 
dynamic to derive relationships between large-scale 
climate variables and the surface climate at point 
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Figure 2. Start date variation during I992-2004 at the 14 ,tudy localities 
(day of the year). Linear trend equation in a \-v represented. The equation 
values and R 2 coefficient for each locality are also indicated. 

locations [47]. The climate scenario applied for the end of 
the 21st century was the x2 C0 2 proposed by the 1PCC. 
which assumes a 1% per year increase in C0 2 concen- 
trations, which is the most probable expected increase for 



the next years, tn this scenario, the HadCM3A2a Regional 
Circulation Model (ROM) developed by the Hadley 
Meteorological Centre provided us with future climate 
data. The RCM provided a horizontal resolution of 50 x 
50 km grid cells, which allowed us to obtain a different 
climatic data base at all study sites (Fig. 1). 

RESULTS 

Pollen season start and Climate Change. The pollen 
season start date variation from 1992-2004 is shown in 
Figure 2. In this figure, the trend over the years is also 
plotted for each site. At all sites, the regression trend 
equation presented a negative slope coefficient, indicating 
that the start of the flowering season tended to be earlier 
in later years. This characteristic was more evident in 
inland central sites such as Madrid and Jaen. Due to the 
considerable year-on-year variation in start dates, equa- 
tions are less well adjusted to a linear model, R 2 coeffi- 
cients not being as close to 1. There was an evident 
gradual advance of the pollen season, which was more 
obvious in inland than in coastal areas, due to a greater 
warming in these study areas. Analysis of climate 
databases (National Institute of Meteorology) shows that 
temperatures recorded in the Iberian Peninsula have 
increased considerably. Figure 3 shows average tempera- 
ture registered in the pre-flowering season (January- 
March), from 1951 to the present, at all study sites. The 
positive slope of the trend equations confirmed the 
increase in temperature. 

Using the GDD° model [26, 27.1, it is possible to assess 
starting dates for years where climate data, but no pollen 
observations, are available. To ascertain the potential 
impact of future climate change on Mediterranean forest, 
future climate data, HadCM3A2a RCM, were applied to 
the GDD° model for the years 2025, 2050, 2075 and 
2099. The predicted increase in temperature- in the Iberian 
Peninsula varies among sites (Fig. 1). These differences 
are likely to give rise to marked geographical variations in 
the advance of phenology (Tab. 2). The predicted mini- 
mum advance will be in northwestern sites, such as 
Santiago, while in inland cities like Leon and Madrid 
flowering could take place up to 1 month earlier. In Table 
2 it is possible to compare the detected advance from the 
beginning of the 90' s until today, with the expected 
advance for the end of the present century. 

Daily pollen concentrations. The curves showing 
daily pollen counts at the various study sites are shown in 
Figure 4. Generally speaking, the Quercus pollen season 
in the Iberian Peninsula lasted 3-4 months, depending on 
the vegetation and climate characteristics of the years and 
areas studied. Several peaks were observed in some cases, 
mostly due to the consecutive flowering of different 
species (see section: Description of Sites). Nevertheless, 
external factors such as spring rainfall or changes in wind 
direction or speed may also have contributed to a more 
irregular curve. 
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. columns. 




Site Begi 


nnin 90's Be innii 
date (d.y.) start dal 


lg 21st Be; 


tinning 21st 
[inning 90' s 
(days of 
difference) 


forecasts!^ 
date (d.y.) 


forecast^ foreca 
date (d.y.) daft 


2075 

'(d.y.) 


,099 

forecasts beginnmj 
differ. 


2009 
cm' 'I 


Santiago 


82 


80 


_ 2 


80 


79 


77 


77 




Vigo 


74 


65 


-9 


55 


50 


36 


32 


-33 


Ourense 


73 


69 




54 


48 


33 




-40 


Barcelona 


90 


80 


-10 


77 


76 


74 


75 




Girona 




85 




80 


75 


70 


72 


-13 


Lleida 




84 




80 


77 


75 


71 


-13 


Tarragona 




90 




80 


76 


76 


75 


-15 




115 


105 


-10 


90 


64 


47 


31 


-74 


Madrid 


106 


82 


-24 


65 


47 


27 


26 


-56 


Cordoba 


90 , 


80 


-10 


71 


58 


31 


23 


-57 


Priego 


91 


75 


-16 


59 


50 


25 


25 


-50 


Jaen 




71 




50 


49 


27 


28 


-43 


Granada 


82 


75 


-7 


64 


54 


31 


30 


-45 


Malaga 


81 


73 




68 


54 . 


43 


34 






88 


79 


-9 








45 


-34 



Similarities were detected between curves for sites 
located in the same climate region, probably due to a 
similar species spectrum and similar meteorological 
conditions. It was therefore possible to distinguish 4 
general patterns in the NW, NE, Central, and Southern 
areas. Sites located in NW Spain (Figs 4a-4d), displayed 
low pollen concentrations with 1 peak right at the start of 
the season and another at mid-season. In the NE area, a 
lower peak is detected at the beginning of the season and 
another higher peak at the end (Figs 4e-4h). Very 
different pollen counts were recorded in the Central sites 
of Madrid (Fig. 4i) and Leon (Fig. 4d), Madrid values 
being substantially higher. Nevertheless, a similar pattern 
with 2 peaks (the higher one at the end of the season) was 
evident. Figures 4j-4n, showing Quercus pollen curves in 
Southern Spain, unhealed high pollen counts, with several 
peaks during the season. 

Daily pollen forecasting. Spearman's correlation test 
detailed results are not shown, due to their length and 
similarity to regression results. They revealed that 
temperature was the main factor affecting daily pollen 
concentrations at all sites. The correlation with 
temperature was positive in the Pre-Peak period, but 
negative in the whole-season analysis. These results were 
particularly evident at sites with longer pollen seasons. 
Moreover, the correlation analysts of pollen vs. wind data 
showed significant positive correlation coefficients at 
coastal sites on days when the wind came from inland; in 
contrast, the wind coming from the ocean prompted a 
decrease in airborne pollen. Similar!}, in inland areas 
where Quercus pollen is influenced by the effect of 
valley-mountain winds, or where Quercus populations are 



located quite far from the monitoring sile, wind direction 
and intensity also influenced pollen counts. Regression 
equations for many sites showed that daily temperature 
also exerted a positive effect on pollen data in the Pre- 
Peak period, while the results for the After-Peak period 
show that humidity and rainfall has a negative effect on 
airborne pollen. Moreover, previous pollen counts are 
more effective forecasting variables in the After-Peak 
period (Tabs 3, 4). Fewer year-on-year differences were 
observed at coastal sites than in inland regions. 

Validation of daily pollen forecasts. Most of the 
proposed daily pollen forecasting models recorded signi- 
ficant statistical coefficients i labs 3. 4). Nevertheless, 
models were externally validated for 2001 and 2004, both 
graphically and statistically. Figure 5 shows the expected 
vs. observed daily pollen values in 2001 and 2004 for 
each site. Grey lines represent actual data and black lines 
represent expected data. VVulei line represents 2001 data. 
The R 2 coefficient and the Wilcoxon matched pair test 
results from the statistical validation are shown in Table 
5, p values > 0.01 indicate non-significant differences. 
The acceptable statistical results and similar patterns of 
the curves suggested that main daily pollen characteris- 
tics, such as the increase, peak or decrease periods, can be 
forecasted using the proposed equations. 

Pollen Index forecasting. Annual Pollen Index (PI) at 
most sites were very high compared to pollen emitted by 
most of the anemophilous species presented in the study 
sites, particularly Pis detected in the Mediterranean area 
(Fig. 6). The Atlantic sites located in NW Spain had the 
lowest Quercus PI; higher amounts were recorded at 



Quercus pollen season response to Climate Change 



217 



Table 3. Daily pollen forecasting equations fr 



lies that suffered the drought period (1992-1995). 



BARCELONA 94-95 


Pre R 2 0.36 p<0.00 
After R 2 0.49 p<0.00 


LEON 94-95 


Pre R 2 0.82 p<0.00 
After R 2 0.53 p<0 00 


MADRID 93-95 


Pre R 2 0.56 p<0.00 
After R 2 0.75 p<0.00 


CORDOBA 93-95 


Pre R 2 0.68 p<0.00 
After R 2 0.60 p< 0.00 


PRIEGO 94-95 


Pre R 2 0.89 p<0.00 
After R 2 0.70 p<0.00 


GRANADA 92-95 


Pre R 2 0.75 p<0.00 
After R 2 0.73 p<0.00 


MALAGA 92-95 


Pre R 2 0.70 p<0.00 
After R 2 0.77 p<0.00 


JALN 96-00 


Pre R 2 0.65 p<0.0() 


BARCELONA 96-00 


Pre R 2 0.57 p<0.00 
After R ; 0.52 p<0.00 


LEON 96-00 


Pre R 2 0.67 p<0.00 
After R ; 0.82 p<0.00 


MADRID 96-00 


After R- 0.49 p<0.00 


CORDOBA 96-00 


Pre R : 0.75 p<0.00 
After R ! 0.59 p<0.00 


PRIEGO 96-00 


Pre R 2 0.94 p<0.00 
After R 2 0.83 p<0.0() 


GRANADA 96-00 


After R 2 0.65 p<0.00 


MALAGA 96-00 


Pre R : 0.55 p<0.00 



QUER: -9.9+0.17Q2+0.98Tmx 

QUER: 20.08+0.1 lQ2-0.94Rf-0.83Rf-0.83Tmn 

QUER: 26.76+7.19WF3+3.85Rfl+2.17WFl 

QUER: 28.3+0.27Ql-0.29H+0.32Rfl -0.24Tmx2 

QUER: -61.1 8+6.33Tmn2+0.3 1 Q1+U7S5 

QUER: 1 3.5+0.4 1Q2-0.23Q1-0. 16Q3 +0.43H3-0.58Tmn2 

QUER: -61.1 8+6.33Tmn2+0.3 1Q1 + 1.17S5 

QUER-108.19+0.68Q1+1.6WF3-1.16WF4+2.18S2 

QUER: -23 1 .2+2.05Q 1 +46.7WF4+4. 1 S3 

QUER: 71 .2+0.47Q1 +0.14Q3-0.52H- 1 .5Tmn2 

QUER: -14.4+0.54Ql+L27Tmxl+0.1 1Q2 

QUER: 40.9+0. 1 9Q2+0.29Q1 -0.48H-0.84Tmn 

QUER: -192.7+0.178Q3+6.21Tmed+0.1 1Q2 

QUER: 4O.9+0.2Q2+0.3Ql+1.6WV-3.48Tmx +0.34WF4 

QUER: -70.9+ 1. 1 4Q 1 +20.4WF1 +9.55WF4 - 1 6.78Rf 

QUER:-17.9+0.24Q2+3.84Tmn2-1.63Rf 

QUER: 254.19+0.98Q1-0.73Q3+0.68Q2 -2.8Tm-22.8S 

QUER: -48.65+0.1 8Q3+3.14Tmx+0.09Q2 -0.13WV4 

.25Q1-0.43WV +0.66H-0.33Tmx2 
I.59Q2-0.6WF4+3.9S1 
D.08Q3-4..5Tm2+0.03Q2 



QUER: 7.98+0. 15Q3-I 
QUER: -4.8+0.51Q1 + 
QUER: 99.2+0.25Q1H 



QUER: 4 
QUER: 4 



,4Hr> 



' -S.4W1-2 



After R 2 0.41 p<0.00 
After R 2 0.55 p<0.00 



()1 IIR: "0 5+ll.7:QI-0.15Q2-6.441mnl 

QUER:-148.09+2.32Q2+0.82WC3+3.33Tmx2 

QUER: 75.39+0.39Ql-0.9Tmn-().79H -0.74Tmn2 

QUER: -305.8+0.6Q2+19.1Tmx-0.2Q3 +)5.8Tmxl+4.2WFl 

QUER: 77.4+0.25Q2-3.05Tmn2-0.3WF2 -2.23Rf 

QUER: 25.6+0.42Q1-0.19H-3.4H1 



Pre - Pre-Peak Season, Quer - Quercus pollen , Tmx - Maximum Temperature, Tm - Mean Tempt 



S - Sunshine Hours, Rf - Rainfall (mm), After - After-Peal 



showed 



, in Medi- 



Velocity (km/h), C - Wind Cal) 
Pollen grains by cubic meter. 

Central and Southern sites. Regression analys 
that rainfall and temperature (particularly 
temperature) were the most important variable: 
terranean areas, while both maximum and uiimiiiuiu 
temperatures prior to flowering were the most important 
parameter in Eurosiberian areas. There was a geogra- 
phical gradient with respect to the months with the 
greatest influence in the forecast. Meteorological condi- 
tions 1 month prior to the pollen season had a decisive 
effect on models (Tab. 7). 

Validation of Pollen Index forecasting. All regression 
models showed an adjusted regression coefficient (R' ) 
close to 1. with a probability value from 95-1001 (Tab. 
7). In Table 6. the plotting of actual vs. expected data 
from the validation in 2001 and 2004 yielded good 
results, with a Root Mean Square Error between actual 
and expected values of 16.2 in 2001 and 21.4 in 2004. An 
acceptable error under 25% was generally found. The 
difference was considerably smaller for Vigo. Ourense, 
Madrid. Cordoba. Granada and Malaga. 



day, Index 2, 3, . 



5, WFn - Wind Frequency, n quadranl. linin 
■s pollen , H - Relative Humidity, WV - Wind 
i from the two, three, ... previous days, gr/m' 



Climate Change model simulation. The average 4°- 
5°C rise in temperature produced by expected global 
warming prompted an increase in modelled flowering. 
Although the model did not predict significant changes in 
rainfall in the Iberian Peninsula, a increase of 5-10% in 
annual rainfall may occur in the Mediterranean area. This 
increase would be registered by means of more torrential 
rainfall episodes alternating with higher drought periods. 
The rate of future climate change simulated by the RCM 
varied significantly among Iberian areas (Fig. 1). Expec- 
ted changes in temperature and rainfall could prompt a 
marked geographical variation in pollen production for 
Quercus species. The use of meteorological parameters 
for 2015. 2050. 2075 and 2090 m the validated Pollen 
Index forecasting models indicated that Quercus pollen 
production in the Iberian Peninsula could become signi - 
ficantly higher at most study sites (Tab. 8). The range of 
variation in PI as a function of the mean PI value (Pollen 
Index/Mean Pollen Index) is shown for the last decade. 
Quercus pollen production would record the maximum 
increase between 2050-2075, while few differences are 



in response to Climate Change 



JlLki 

8 




Figure 6. livoliilion of annual Qmraix Pollen Index in study areas. 

expected for the last 25 years of the 2 1st century. Compa- 
rison with the variation factor expected in the future 
(Future Pollen Index/Mean Pollen Index) shows that 
expected PI will increase two and a half- fold from today's 
values. Climate simulation for the doubled C0 2 scenario 
produced different rates of pollen increase in different 
areas; these were lower in coastal areas than at inland 
sites, Mediterranean inland areas being the most affected. 

DISCUSSION 

Understanding pollen emission dynamics is fundamen- 
tal for the characterization of potential allergens that may 
be of greater health relevance in both natural spaces and 
populated areas. Quercus airborne pollen counts at a 
given sampling station reflect reproductive phenology 
over an extended area [301. Taking into account the 
number and distribution of the sites studied in the present 
paper, this research provides valuable information on the 



pollination dynamics features of Quercus genus in the 
Iberian Peninsula. The present results could be very 
useful for allergy physicians to know the main meteorolo- 
gical factors influencing on Quercus pollination dyna- 
mics. Moreover, the investigations carried out on the 
possible future consequences of climate change can give 
them a new vision of the possible changes in atmospheric 
allergens occurrence. 

Results indicate that Quercus pollen season start has 
been advancing over recent years in the Iberian Peninsula. 
Results were similar to those obtained in a previous study 
involving fewer sites and years [27]. Other authors in 
northern Europe have reported similar findings for the 
Quercus pollen season [9]. With regard to other species, a 
progressive advance in flowering has also been observed 
in Betula [17J, Poaceae [18], and Alnus and Olea [11] 
over the last 20 years. All these findings indicate that 
spring phenology in Europe is being significantly affected 
by global warming. On the basis of previous studies [26, 



Table 4. Daily pollen forecasting equation-; ... 
from localities with a shorter data base, iGiron; 



orth-west localities that did not surfer the drought period (Santiago. Vigo and Ourense). Results . 
i and Tarragona i also shown. 



SANTIAGO 93-00 


Pre R 2 0. 


4 p<0.00 


QUER: 


26.9+0.56Q1 


e(>. 1Q3 + 


.9Tm+1.8Tm.\l 




After R- 




QUER: 


16.7+0.5 1Q1- 


1.41H-0.3 


2C+0.51Tmx2 


VIGO 95-00 


Prc R" 0. 


M p<0.00 


QUER: 


7.69+0.38Q2 


-n.UiQI 


).9WF3-1.59S2+2.7Tmx2 


After R 2 


1.68 p<0.00 


QUER: 


7.61+0.28Q1- 


2.1WE1 -11.1150* l.dOWId ti .MRU 


OURENSE 93-00 


Pre R : 0. 




QUER: 


26.9+0.72Q1 


s-1.5Tm,\ 


+-0.9S 




After R 2 


).81 p<0.00 


QUER: 


J.8+0.50QM 


O9H-0.3 




GIRONA 96-00 


Pre R 2 0 




QUER: 


-168.1 1 + I.07QI + I5.3 1 






After R 2 


).73 p<0.00 


QUER: 


94.2+0.71Q1 


0.14Q2- 


..5Tmn2-1.2S2+1.08WV 


LLEIDA 96-00 


Pre R 2 0 


54 p<0.()0 


QUER: 


-207.5+O.89Q 


1+9.68T 






After R 2 


0.51 p<0.00 


QUER: 


12.49+0.58Q 


-3.6Tmn 


l-0.08Tmx+1.72Rf3 


TARRAGONA 96-00 


Pre R 2 0 


53 p<0.()0 


QUER: 


-44.04+6. 46T 








After R 2 


0.52 p<0.00 


QUER: 


I9.86+0.59Q 


+0.59Rt 


-0.35H+0.19H1 
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Site 




R 


Wilcoxon Mat 


die.! Pair.- 


. lest 




N 


T 


z 


plevd 


Saul uigo 
















2004 


068 


122 




1 24 


0 21 


Vigp 








1 ~J>(! on 


1 42 


0 15 




"-II0 4 


046 


107 




7 85 


000 


Ourense 






116 










-IK 11 


065 




'"'ll'j -° 


077 


044 


Leon 








^900 








2004 


072 


81 






000 


Barcelona 












On' 




2004 


064 


106 


o'oor 


1 83 




( Tin ma 


^004 


0.73 


121 


8 s.. .SO 


3.65 






J II 14 






LscUcOO 




0 si 


Tarragona 






'99 










HH)t 


083 




54550 


.23 


026 


Lleida 






122 


l.-SmOO 








2004 


061 


106 


Lc.37.s0 


■1 68 


009 


Madrid 






12^ 


1,201.00 




I! 04 




2004 


034 




2,224.50 


038 




Cordoba 


2001 


078 






1 21 






2004 


0.65 


98 


1.306.00 


2.33 


0.0s 




2001 


0.58 


112 


1.309.00 


5.00 


0.01 




2002 


0.45 


77 


LOW. 50 


1.55 


o.i? 




2001 


0.79 


91 


783.00 


1.16 


0.2 1 




2004 


0.23 


104 


184.00 


4.14 




( jrnnada 


2001 


0.79 


127 


1.033.00 


0.05 


0.96 




2004 


0.43 


119 


493.00 


1.97 


0.06 


Malaga 


2001 


0.52 


131 


1,407.00 


1.53 


0.13 




2004 


0.40 


1 15 


2.839.00 


0.64 


0.52 



Table 6. Annual Pollen Index 1 
The Expected (Exp) and Actual 
between them are shown. The fi 
error (%Error) in relation with 



validation for 2001 and 2004. 
)llen Index (PI), and difference 
mn indicates the percentage of 
al Pollen Index. The Absolute 



1 Marked differences are signilm am al p. 0 (is 

27] it is clear that the Quercus start date in the Iberian 
Peninsula is strongly determined by temperature, there- 
fore, if our climate continues to change as projected, the 
start of the Quercus pollen season would continue to 
advance. 

Comparative aerohiological studies of sites with diffe- 
rent climatic characteristics give us a global idea of the 
different behaviour of airborne pollen under specific 
meteorological conditions 111 large areas I 10). Regarding 
Quercus daily pollen dynamic in Spam, the seasonal 
floral phenology of the various Quercus species presented 
at each she determines the shape of the pollen curve. The 
curve was fairh similar for the ^ stuck sites in the NW 
area: Santiago, Vigo and Ourense. In Southern Spain, the 
3 species present in lover altitude areas 1 Cordoba and 
Malaga). Q. ilex subsp. ballota. Q. cocci/era and Q. suber 
have consecuthe flowering periods that gave rise to a 3- 
peak pattern, the highest being the Q. ilex subsp. ballota 
peak. In higher altitude areas, such as Jaen. Granada and 
Priego, the great abundance of Q. ilex subsp. ballota 
accounted for almost all pollen grains reflected in the 
citne. In these areas, Q. ilex subsp. ballota populations 



Site 


Year 


Expected 
PI 


Actual PI 


Dif 

(Exp- Act) 


1 11... 


Santiago 


2001 


209 


282 


-73 


25 




2004 


792 


685 


107 


15 


Vigo 


2001 


520 


629 


-109 


17 




2004 


1,972 


1.742 


-230 


1 


Ourense 


2001 
2004 


1,386 
2,168 


1,387 
2.578 


410 


0.07 
15 


Barcelona 


2001 


2,424 


5,965 


-3,541 


59 




2004 


4,289 


6,002 


-1,713 


28 


1 arragona 


2001 


1,410 


1.412 


-2 


0.14 




2004 


2,647 


3,931 


-1,284 


32 


Girona 


2001 


1,837 


1,411 


426 


30 




2004 


20,933 


19.113 


-1,820 


9 


Lleida 


2001 


1.122 


LI 19 


3 


0.26 




2004 


3.225 


5,429 


-2,204 


40 


Leon 


2000 


1,197 


822 


-375 


41 




2004 


1,244 


1,405 


-161 


11 


Madrid 


2001 


12.564 


10,307 


2,257 


19 




2004 


2,422 


2,525 


103 




Cordoba 


2001 


13,324 


14.260 


-936 


6.5 




2004 


2,990 


2,891 


-99 




Priego 


2001 


8,024 


7,774 


-250 


3 




2002 


4,272 


6,158 


-1,886 


30 




2001 


5,931 


7,161 


-1,230 


17 




2004 










Granada 


2001 


5,235 


4,429 


806 


17 




2004 


3,992 


4,303 


-311 


7 




2001 


6,224 


8,156 




23 




2004 


3,195 


2,823 


372 


13 




2001 
2004 






852 
838 


17.4 
18.4 


SD 


2001 
2004 
2001 
2004 


4,367 
5,026 


4,372 
4.565 


16.2 
21.4 





t different altitudes caused consecutive peak 
1 patterns: in this case, topography was respon- 



pollin 

sible for the irregular curve. 

The specific aerohiological behaviour of Quercus 
pollen has been studied in different parts of Europe and 
the USA from a descriptive standpoint [9, 25, 35, 37, 501. 
However, no studies have been performed on the 
forecasting of airborne Quercus pollen levels. The daily 
forecasting results presented here indicate that the main 
variables were temperature and hours of sunshine, with a 
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Table 7. Annual Pollen Index (PI) forecasting regression equations. 






Site 


. Regression equation 


R ; 


p level 


Santiago 


PI=-22885+ 1 9(RfFl )+460(MeanF2)+ 1 086(MxJ2)+270(MxJ 1) 


0.99 


0.02 


Vigo 


PI= 1 1 345+4(RfJ2)+ 161 0(MnF2 )t-434(MnJ2)- 1652(MxF2) 


0.99 


0.02 


Ourense 


PI=2905+5.4(RfJl)+3.6(RfJ2)+(MnFl)-274(MnJl)-238(MeanF2) 


0.99 


0.03 


Barcelona 


Pl=440 1-61 (R IF 1 )+28(RfJ2)-l 50(Rf2) 


0.99 


0.00 


Tarragona 


PI=36344-2423(MxE2)+228(RfFJ) 


0.99 


0.00 


Girona 


PI=3522 1 -2023(MxM 1 )+ 1 238(MeanF2)-76(RfJ 1 )+4()(RfJ2) 


0.99 


0.00 


Lleida 


PI=5075 -45 ( Mx M I )- 1 39(MnM 1 )- 1 2 1 (RIM 1 )+205(RfFl ) 


0.99 


0.00 




PI=9047-888(MxJ2)+8(MxF2)+5.3(RF 1 ) 


0.98 


0.00 


Madrid 


PI=-26098+ 1 955(MnF 1 )+3909(MnM 1 )+996(MeanJ 1 )+1005(MeanJ2) 


0.99 


0.02 


Cordoba 


PI=32909+ 1 652(MxF2)+ 1082(MeanF 1 )+ 1 00(RfFl)+4 1 5(MnE2) 


0.98 


0.00 




PI=757+1367(MnF2)+52(RfF2) 


0.99 


0.01 




PI=- 37045-357 1 (MeanM 1 )+6664(MnM 1 )+239 1 (MxJ 1 ) 


0.99 


0.02 


Granada 


PI=3731-1002(MeanFl)+1002(MnFl)-64(RfMl)-21(MeanF2)+10.6(RfJl)-170(MeanJl) 


0.99 


0.05 


Malaga 


PI=5576+648(MxF2)+9.4(RfJ 1 )+409(MeanF2 ) 1 303( Mx M 1 I 


0.99 


0.00 



RfJ - Rainfall January , RfF - Rainfall February , RfM - Rainfall March , MxJ - Maximum Temperature January, MxF Maximum Temperature 
February, MxM - Maximum Temperature March, index 1 1st fortnight. MeanJ - Mean Temperature January. McanF - Mean Temperature February. 
MeanM - Mean Temperature March. Mn.i - Minimum Temperature January. MnF - Minimum Temperature February, MnM - Minimum Temperature 
March, Index 2 - 2nd fortnight. 

Table 8. Average Pollen Index (PI) from last years (1992-2004) compared with the calculated future PI for 2099 at different sites of the Iberian 
Peninsula. Range and average value of variation factors from the last decade is compared with the adjusted increase factor expected lor the end ot the 
21st century. 



Site 


Last decade 
Mean Pollen 
Index (MPI) 


Last decade 
Range of 
Vai ration 
Factor 

(PIln/MPI) 


Last decade 
Average 
Variation 
Factor 


PI 2025 


PI 2050 


PI 2075 


PI 2099 


Variation 
Factor 
(PI 2099/MPI) 


Santiago 


825 


(0.3-2.5) 




870 


870 


880 






Vigo 


2,044 


(0.4-2.2) 


1.3 


2.100 


2,100 


2,200 


2,250 




Ourense 


1,500 


(0.9-2.2) 


1.5 


2,050 


5,600 


9,650 


10,800 




Barcelona 


4,500 


(0.4-1.2) 


0.8 


5.100 


5,900 


6,200 


6,100 




Tarragona 


3,000 


(0.7-1.4) 


1.05 


3,510 


4,100 


4,650 


4,640 


1.5 


Girona 


9,850 


(0.2-1.7) 


0.95 


10.530 


13,730 


14,800 


15,410 


1.5 


Lleida 


3,500 


(0.4-1.7) 


1.05 


5.080 


7,550 


7,900 


8,000 


2.3 


Leon 


900 


(0.5-2.3) 




2,500 


3.740 


3,970 


4.100 


4.5 


Madrid 


7,500 


(0.4-2.0) 


1.2 


9,000 


16.750 


21,200 


23,800 




Cordoba 


8,000 


(0.3-2.3) 




11,210 


16.600 




24.500 






6,000 


(0.3-2.4) 


1.35 


6,400 


10,500 


17,320 


17,600 






4,400 


(0.3-1.8) 


1.05 


6,500 


7,540 


8,500 


9,800 


2.2 


Granada 


3,200 


(0.5-1.3) 


0.9 


4,870 


5,600 


5,600 


5,600 




Malaga 


5,000 


(0.4-1.6) 


1 


5.090 


5,500 


6,210 


6,600 




Absolute Mean 


4,300 


(0.4-1.9) 


1.16 


5,340 


7,580 


9,470 


10,005 





MIT - Mean Pollen Index. PI Pollen Index. IT 1 n - Pollen Index from each study year. 



positive effect on daily Quercus pollen variation. It was 
noticeable that at many sites, daily temperature variables 
showed positive coefficients in the Pre-Peak equations, 
whereas previous results for the whole season displayed 
significant negative coefficients. The extensive length of 
the Quercus pollen season makes coincident the increase 
in temperature in late spring/early summer with the end of 
the season, when decreasing pollen concentrations were 



recorded due to the intrinsic phenology of the plants. Due 
to this fact, different meteorological factors influenced in 
the 2 pollination periods. Results showed that different 
meteorological factors affected in the 2 periods. In the 
Pre-Peak period, temperature variables positively affect 
the pollen liberation, whereas in the After-Peak period. 
Humidity and Rainfall influence the atmospheric dyna- 
mics of the pollen already liberated. 
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Moreover, the use of previous pollen records has 
proved to be a very useful tool in pollen forecasting 
overall in the After- Peak period. The use of this sort of 
data is recommended in taxa with a recurrent pollen curve 
year after year, such as Quercus [6, 43]. In general, the 
validation of the proposed models was performed in 2 
years in order to confirm the results. The \ears from 
2001-2004 were climatically similar to the 1996-2000 
period, and drought years were not detected. 2001 and 
2004 years were randomly chosen. Validation of the daily 
forecast models lor 2001 and 2004 showed that the main 
daily pollen characteristics can be accurately forecasted: 

A number of authors have studied inter-annual Pollen 
Index variations in different species in broad geographic- 
zones, such as all Europe [55], England [18], Italy [21| 
and the Iberian Peninsula [25, 34]. In general, these 
authors reported that bio-geographic and climatic 
characteristics, as well as hun'lan action, can affect air- 
borne pollen counts in each area. In the present study, the 
highest pollen concentrations were detected in different 
years, depending on the area, although fewer inter-annual 
differences in pollen counts were observed at coastal sites 
than in inland study areas; this may have been due to the 
lower degree of climate variation in these areas. In 
genera], the annual PI tended to decrease in the early 
years of sampling (1992-1995) due to the drought. Later, 
annual PI increased due to recovery of these trees. 
Nevertheless, this recovery was not detected in the first 
"post-drought" year (1996) as might be expected, but 2 
years later, in 1997. The extreme drought period 
weakened oaks, slowed down vegetative growth, and 
prompted biomass loss. The wetter 1995-1996 winter 
probably led to a recovery of vegetative biomass and to 
strong production of vegetative buds instead of floral 
structures. This phenomenon was clearer in the south of 
the Peninsula, where the shortage of rain was most acute. 
The same trend has been observed in other species in the 
Peninsula, such as Olea europaea [13| or Platanus spp. 
[14]. 

In 1997, 1998, 2003 and 2004 the highest PI was 
recorded in the southern cities, whereas in the northwest 
area - unaffected by the drought - the highest 
concentrations were detected in 1995. During this year, a 
very high concentration of Quercus pollen was also 
observed in England [9j. Weather conditions over 
previous months have been shown to exert a remarkable 
effect on the intensity of Quercus flowering. Rainfall 
during the month prior to the pollen season has been 
highlighted as the parameter making the greatest 
contribution to pollen production. In other species, such 
as the olive which flowers in late spring, studies report 
that in Cordoba (southern Spain) rains make the greatest 
contribution to flower production in March [23]. whereas 
in Perugia (Central Italy) the ram effect on flowering 
intensity is strongest in May |20|. In the Cupressaceae 
family which flowers in winter, it has been reported that 
the minimum temperature in November and rainfall in 
November and December were the variables everting the 



greatest influence on the PI [22]. In other models 
designed for trees with clear biennial patterns, such as 
Betala, PI values from the previous year have been 
included as independent v ariables, as well as temperatures 
in the current year [12]. In the results obtained here, the 
variables with greatest influence on total Quercus pollen 
counts were rainfall and maximum temperature in 
Mediterranean localities, and minimum temperatures in 
Eurosiberian Region sites or these influenced by its 
climate. This may be due to the fact that Mediterranean 
species flower during early spring when minimum 
temperatures do not vary significantly from winter, 
whereas Eurosiberian species flower in late spring when 
minimum temperatures are much higher and go a long 
wav towards determining the development of floral 
structures. The gradients for the most influential months, 
i.e. weather characteristics 1 month prior to flowering, 
have the strongest influence on oak flowering intensity. 

In Spain. 21.5% of the population suffer from pollen 
allergy which involves a high number of sick leaves [2]. 
Annual allergy costs for the Spanish Government rise 
above 1.500 million Euros |39|. Advance knowledge on 
the intensity of the pollen season would help allergy 
physicians and allergy sufferers to better plan and manage 

Prentice et al. [49] examined the implications for future 

concentration by 2099 ranged from 1.9-3.5 times. Never- 
theless, the most expected C0 2 increase is projected to 
rise 2 times [10]. The implementation of future climate 
data presented in this work was carried out under x2 C0 2 
scenario. Results indicate that the Quercus species in the 
Iberian Peninsula are already responding to climate 
warming. Changes in temporal occurrence of atmospheric 
allergens could have important clinical consequences [3]. 
The floral phenology of Quercus species in the lberian 
Peninsula is suffering the rising spring temperatures in the 
Mediterranean area, like the phenology of Betula and 
Almts in colder biomes, where climatic warming is 
already having a visible impact on vegetation [8, 16, 41, 
54]. Moreover, floral phenology of other tree species in 
the Iberian Peninsula is being affected. It has been 
detected that Olea flowering season has suffered a general 
advance in the Iberian Peninsula [24|. This advance has 
been shorter than that observed for Quercus species, 
probably due to the higher dependency on temperature of 
early spring tree species flowering date. For late spring 
species, such as olive, photoperiod and rainfall also 
influence the flowering date [24]. The increase of 
temperature is probably affecting the flowering of many 
other Iberian species, especially species like Platanus or 
Ulmus. which have shown a clear response to temperature 
before flowering |1]. Modelling future Ombrothermal 
Indexes combining Temperature. Rainfall, Aridity and 
Continental Indexes for the inland central area of Leon 
suggest that the next 50 years will not see significant 
changes in seasonal rainfall distribution - a major factor 
for climate-related plant species [51]. The study also 
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predicts an increase in temperature of 1-2°C over the next 
50 years, which will lead to accelerated vegetative and 
flower growth and an advance in flowering dates. The 
results of that local study agree with those reported here 
for the Leon area. 

A number of authors have reported that the increase in 
airborne pollen counts recorded by most European 
aerobiology stations over recent years ma\ be due to this 
climate warming [33, 47]. Some authors, using a 
vegetation biogeochemical model (BIOME3) where 
torrential rainfall increases 5-10',-r in the Mediterranean 
area, indicated that future climate conditions will lead to 
the spread of an evergreen forest in the Mediterranean [7, 
31]. This hypothesis, taken in conjunction with the 
models proposed here, suggests an increase in the 
evergreen Quercus pollen of the Iberian Peninsula. The 
pollen season of these trees will start earlier, and will also 
prompt higher pollen emissions throughout the Mediterra- 
nean area. These effects would be more evident in inland 
areas as Madrid, Leon, Jaen or Cordoba. 

The present paper, together with other studies 
performed in different anemophilous species, show 
considerable evidence to suggest that climate change, 
specially temperature increase, will have, and has already 
had, impacts on aeroallergens. Climate change effects 
include impacts not only on pollen amount and pollen 
season, but also in plant and pollen distribution and other 
plant attributes. For this reason, theoretical models, such 
as these here presented, must be understood only as a 
partial approach to the future environmental conditions 
that could also be modified by deforestation, or land use 
changes. There are many research challenges to a more 
complete understanding of the impacts of climate change 
on aeroallergens and allergic diseases [3]. In-depth studies 
covering main allergenic pollen species in Europe should 
be carried out in order to obtain a full idea of the possible 
effect of climate change on the pollination dynamic 
characteristics. Therefore, new efforts on this subject 
should be made by public health authorities and allergy 
physicians in order to be aware of these changes in the 
environment to adapt allergy treatments. 

CONCLUSIONS 

Quercus floral phenology has been advancing over the 
last few years; pollen seasons in the Iberian Peninsula 
have, since the beginning of the 90s, been starting an 
average of 2 weeks earlier. 

Quercus pollen-curve characteristics varied between 
different bio-geographical /ones, but were similar for 
sites with similar climate and species. It was therefore 
possible to distinguish 4 general patterns in the NW, NE, 
Central, and Southern areas. In the Mediterranean area, 
minimum and maximum temperatures were the main 
parameters influencing daily pollen-count. Maximum 
temperature was the most influential factor in the 
Euroiberian area. 

The Pollen Index varied considerably between sites and 
years, depending on previous weather conditions. In the 



Mediterranean area, rainfall and temperature (particularly 
maximum temperature) prior to flowering were the most 
important variables. Both maximum and minimum tempe- 
ratures were the major parameters in the Euroiberian area. 

Results of external validation of daily and annual 
forecasting models suggest the possibility to use Quercus 
pollen forecasts in future in the Iberian Peninsula. 

Using the RCM climate-change model, it may be 
suggested that the impact of future climate variation on 
Iberian Quercus floral phenology would be different 
depending on the areas. It would be more * marked in 
inland than in coastal areas. In general, the pollen season 
will start earlier and will also prompt higher pollen 
emissions throughout the inland Mediterranean area in 
general. Tims, the likelihood is that climate change will 
modify the time of the occurrence of human oak-pollen- 
allergy problems. 

Using climate predictions from the RCM climate, our 
results would suggest that, as climate warms, the Iberian 
Quercus pollen season will start earlier and will also 
prompt higher pollen emissions, with the strongest effects 
in inland areas. Therefore, the likelihood is that climate 
change will intensify human oak-pollen-allergy problems. 
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